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NOMENCLATURE 


lateral cyclic pitch measured from hub plane and in **wind-hub** system, 
rad (or deg) 

lateral cyclic pitch measured from hub plane and in hub-body system, rad 
(or deg) 


blade lift-curve slope 

longitudinal first-harmonic flapping coefficient measured from hub plane 
and in ”wind-hub’* system, rad 

longitudinal first-harmonic flapping coefficient measured from hub plane 
and in hub-body system, rad 

blade coning angle measured from hub plane, rad 

longitudinal cyclic pitch measured from hub plane and in "wind-hub** 
system, rad (or deg) 

longitudinal cyclic pitch measured from hub plane c*nd in hub-body system, 
rad (or deg) 


lateral first-harmonic flapping coefficient measured from hub plane and 
in "wind-hub** system, rad 

lateral first-harmonic flapping coefficient measured from hub plane and 
in hub-body system, rad 

H force coefficient, C 


^ PTTR2(fiR)2 

Qx. 


torque coefficient, = 


piTR2(nR)2R 


thrust co( fficient, = 


pirR2(nR)2 

Y force coefficient, C,, = r 

Y pirR2(fiR)2 

blade chord, in 


flapping hinge offset, m 

component of main rotor resultant force in the rotor disc plane in 
i|^* = 0 direction, N 

blade moment of Inertia about flapping hinge, kg-m^ 



Hi 


shaft tilt w.r.t. fuselage, positive forward, deg 


Kl 




k2 




L 


F 




P 

P 

Pw 

• 

P 

^F 

q 

c 

q 

R 


pitch-flap coupling ratio, tan 63 
flapping hinge restraint, N-m/rad 

desired longitudinal control sensitivity, kj 


3a 


Is 


36 , 


desired lateral control sensitivity, k2 




main rotor contribution to aircraft rolling moment about body axis x 
N-m 


>1 


rolling hub moment in hub-body system 
main rotor hub moment about Xg axis, N-m 

main rotor contribution to aircraft pitching moment about body axis y 
N-m 

pitching hub moment in hub-body system 
main rotor hub moment about yg axis, N-m 
blade mass moment about the flapping hinge, kg-m 
number of blades 

main rotor contribution to aircraft yawing moment about body axis 2. 
N-m 


ratio of flapping frequency to rotor system angular velocity 
aircraft roll rate, rad/sec 

aircraft roll rate in wind-hub system, rad/sec, P„ = P cos S„ + q sin 
aircraft roll acceleration, rad/sec^ 
main rotor torque about Zg axis, N-m 
aircraft pitch rate, rad/sec 

aircraft pitch rate in wind-hub system, rad/ sec, q^ - -p sin 6,^ + q 
aircraft pitch acceleration, rad/sec 
rotor radius, m 


Iv 




aircraft yaw rate, rad/sec 

radial station of the blade element measured from the flapping hinge, 
main rotor thrust force acting perpendicular to rotor disc plane, N 
true airspeed, m/sec 




, m/sec 


uniform induced velocity, v. * ;; TTTTi 

^ 2(y2 + 

main rotor force along body axis, x, N 

nondimens ional radial station of the blade element, x 


A e + r’ 


m 


wind-hub system 


hub-body system 

component of main rotor resultant force in the rotor disc plane in 
= 90® direction, N 

main rotor force along body axis z, N 
hub plane angle of attack, deg 

blade flapping angle measured from hub plane, rad (or deg) 

rotor sideslip angle, that is, the angle between Xg and x^, deg 

Lock number, - 

6 

blade mean profile drag coefficient 
lateral control displacement 
longitudinal control displacement 
e/R 

blade pitch angle measured from hub plane, 

0 = Oq - cos ij; - sin + x0^ - K^B, rad (or deg) 


V 



6q blade- root collective pitch measured from hub plane, rad 


0^. total blade twist (tip with respect to root), deg 


P 

0 


r 


inflow ratio. 


V sin a - V. 


advance ratio. 


QR 

V cos a 




air density, kg/m^ 
rotor solidity ratio 
control advance angle 


azimuth angle measured from downwind in the sense of rotor rotation, rad 
(or deg) 


azimuth angle measured from -Xg in the sense of rotor rotation, deg 
rotor system angular velocity, rad/sec 
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A SIMPLIFIED ROTOR SYSTEM MATHEMATICAL MODEL FOR 
PILOTED FLIGHT DYNAMICS SIMULATION 
Robert T. N. Chen 
Ames Research Center 


SUMMARY 


This report documents a simplified analytical mathematical model of the 
helicopter main rotor; the model was developed primarily for real-time pllot- 
in-the-loop Investigation of helicopter flying qualities. The mathematical 
model explicitly includes the tip-path plane dynamics and several primary 
rotor design parameters, such as flapping hinge restraint, flapping hinge 
offset, blade Lock number, and pitch-flap coupling. The model has been used 
in several exploratory studies, recently performed at Ames Research Center, of 
the flying qualities of helicopters with a variety of rotor systems. 


The report describes the basic assumptions used and the major steps 
involved in the development of the set of equations listed. The equations 
consist of the tip-path plane dynamic equation, the equations for the main 
lotor forces and moments, and the equation for control phasing required to 
achieve decoupling in pitch and roll due to cyclic inputs. 


INTRODUCTION 


There is a need for an expanded flying-qualities data base for use in 
developing design criteria for future helicopters. A safe and cost effective 
way to establish such a data base is to conduct exploratory investigations 
using piloted ground simulation techniques, and to then substantiate the 
results in-flight by using variable-stability research helicopters. 


Mathematical models for real-time, piloted, ground simulation valid for 
specific missions and tasks are therefore needed. Unfortunately, several 
existing general purpose helicopter simulation models, such as the "C-81" 

^ef. 1) and "REXOR" (ref. 2), are not suitable for real-time applications. 
These sophisticated blade-element models are far beyond the capability of 
current-generation computers for real-time simulation. Moreover, because they 
are aimed at multidlsclpllnar:. users, they are inefficient and often do not 
provide the insight needed for designing experiments for Investigating general 
flying qualities. 


Simpler flight dynamic models for piloted simulation exist (refs. 3 4) 

specific existing helicopters. Also, these models are quasi- 
static in nature and do not explicitly account for the dynamic effect of rotor 
modes, such as rotor-blade flapping, which can be Important in studies of 

y ng qualities. Therefore, there is a need to develop general but simplified 
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flight dynamic models for exploratory studies of flying qualities and flight- >• 

control systems and, especially, for use on ground simulators that have only 
moderate computational speed and capacity. 

This report presents a simplified general main rotor mathematical model 
that has been used in experimental design analyses and ground-simulator Inves- ' 

tlgations of helicopter flying qualities in visual terrain flying tasks and * \ 

instrument flight tasks at Ames Research Center (refs. 5-7). ^ The model , 

explicitly contains rotor system design parameters, such as flapping hinge 
restraint, effective hinge offset, blade Lock number, and pitch-flap coupling, 
as shown in figure 1. These parameters have been used in appropriate combina- 
tions to cover a wide variety of rotor systems. 



Figure 1.- Main rotor configuration and parameters. 

The report also describes the tip-path plane dynamic equation, the 
development of the main rotor force and moment expressions, and the control- 
mixing equation for decoupling pitching and rolling moments due to cyclic 
control displacement in the cockpit. 


TIP-PATH-PLANE APPROXIMATION FOR FLAPPING DYNAMICS 

Essential to the main rotor mathematical nx)del is the description of the 
blade-flapping dynamics. For the model described in this report, the flapping 
dynamics v'ere approximated using a tip— path plane representation. A detailed 
description of the development was given in reference 8. For present pur- 
poses, a brief summary of the basic assumptions used and the major steps 
Involved In the development Is presented. 


^Mlyajina, K. ; and Chen, R. T. N. : Analytical and Experimental Study of 

an Advanced Stability and Control Augmentation System for a Hingeless Rotor 
Helicopter (in preparation). 
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Basic Assumptions 

The flapping equation of motion of the rotor blade was first developed 
using the following assumptions. The assumptions are similar to those used 
ror the classical" equations (refs. 9 and 10). 

1. Rotor blade was rigid in bending and torsion, and the twist of the 
blade was linear. 

2. The flapping angle and Inflow angle were assumed to be small and the 
analysis utilized a simple strip theory, 

^ effects of the aircraft motion on the blade flapping were limited 

to those due to the angular acceleration p and q, the angular rate p, q. and 
the normal acceleration. 

4. The reversed flow region was Ignored and the compressibility and 
Stall effects disregarded. ^ 


used. 


5. The inflow was assumed to be uniform and no inflow dynamics were 

6. The tip loss factor was assumed to be 1. 


M wi equation of motion for a two-bladed teetering rotor and 

rotors are given in appendix A. They explicitly contain 
. p imary design parameters, such as flapping hinge restraint, hinge offset, 
blade Lock number, and pitch-flap coupling. The blade flapping, 6(t) in 

those equations was then approximated by the first harmonic terms with time- 
varying coefficients; that is, 

7. 6(t) = ag(t) - aj(t)cos i|) - bj(t)sin \p 

The first and second time derivatives of B(t) are therefore, 

S(t) = ag(t) - [aj(t) + b^(t)nlcos ;p - [bj(t) - aj(t)Q]sin 

e(t) = ag - (aj + 2bjfi - a^S -)cos t'' - (b^ _ 2a,i:> - ^ 


Tip-Path Plane Dynamic Equation 

a.(t)"°I U) "and'h'm"® N-bladed rotor, the time-varying coefficients. 
tem^w?th 'a ^ obtained by equating the constant term and the 

I ^ the flapping equation (A2) using the above 

equatior*- ’"^^ult is the following set of tip-path plane dynamic 
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where 


f 

I and 


p2 = 


K, 


eM„ 




1-K, 



') 


T 

5 .” ^ 1 ’ ^ 1 ^ 

For a two-bladed teetering rotor, the tip-path plane representation loses its 
physical meaning. However, if the approximation of item (7) above for the 
blade flapping is employed, then 3 q is to be treated as a preset constant. 

The coefficients (t) and bj(t) can then be solved by setting e = a. = a. = 0 
in equation (1). ^*00 


MAIN ROTOR FORCES AND MOMENTS 



ri 


The main rotor thrust, the H and Y forces, the hub moments, and tue 
torque were developed for the rotor system of interest. In the development of 
these forces and moments, the same set of basic assumptions (1 through 7, 
above) discussed in conjunction with the development of the tip-path plane 
dynamic equations was utilized. Thus, aerodynamically , momentum theory was 
used in conjunction with the uniform inflow; simple strip theory was utilized 
and the blade forces were analytically integrated over the radius. Because 
the reversed flow region and the stall and compressibility effects were 
ignored, the total rotor forces and moments were again analytically obtained 
by summing the contributions, to each blade, that were analytical functions of 
the azimuth. Because of these assumptions and simplifications, the results of 
the analysis are valid only for a limited range of flight conditions. Never- 
theless, a previous study (ref. 10) has shown that this type of analysis is 
valid for stability and control investigations of the rotorcraft up to an 
advance ratio of about 0.3. Also, similar to the development of the tip-path 
plane dynamic equations, these rotor forces and moments were first obtained in 
Che wind-hub coordinate system. They were then transformed into the hub-bodv 
system (see fig. 2). ^ 

The forces and moments thus developed contain periodic terms; the highest 
harmonic terms correspond directly to the number of rotor blades. For example, 
tor a three-bladed rotor, the force and moment equations contain only 3/revolu- 
tion narmonic terms, and for a four-bladed rotor, 4/revolution harmonic terms 
(this IS shown in appendix B). 

With the harmonic terms dropped, the force and moment expressions for a 
nonteetering rotor are as follows. 
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RELATIVE WIND 



DOWNWIND 

x$, Yg. 2g = ROTCR ''HUB-BODY" SYSTEM 
VS' H ^ rotor "HUB-WIND" SYSTEM 

Figure 2.= Coordinate systems for the main rotor. 
Main Rotor Thrust 


The main rotor thrust expression for a nonteetering rotor is: 


I pacR(fiR)2 i (1 - c^)\ + 0 


- I (1 - - K^bi) - ag 


»[i + T <1 - ^ ®t[i + T 

[i + T‘ O - + ai[y e(l - e)| 


(i - e:2) 


^ [U 

\3 2/ Q [4 


(1 - c)' 






Main Rotor H and Y Forces in Hub-Body System 

The expressions for the main rotor H and Y forces in the hub-body 
system are: 


2C 

H„ = § pacR(f2R)2 ~ 
r L ao 
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Main Rotor Hub Moments 


The expressions for the main rotor hub moments are: 

+ W sin 8 ^ 

^HF “ S„ 


where 


<V» ' f [vi - ^ * 2b, i! - a,t!2)l - I 


i ^ T « - 


’’ <*lc - - i (1 - . 

* i(' a ®w + a ‘=‘’® ®v 


H.2^«..,b,.(i-|)(;i.b.) 


N r 1 

2 ^ “ 2 ajfi - bjfi 2 )J _ | - (1 - - K^a^) 

r — 


-[i + I - c>](Blc - >'l'>l) + 3 't + f <1 - E)A t 

X <1 - T * (I - |)(x - =i) ^ i e» + f 

Main Rotor Torque 

The expressions for the main rotor torque are: 

2 c 

Qp = ~ pacR 2 (f^R )2 __S 

^ ^ aa 

^ II * (1 ( 6 , - K.a,)fi Xf - ^ ^ 


y (1 - e)aj 


sin S. 


-K.a,,[X(f-i)TX(T) 

(■a' - *l) I (f »w J 51" B„) - (A,^ - K,a,)j(i - + bj) 

^0 ■*■ i7 ■*■ g(^ Q sin ^ cos 

■c - ^ v[o - f)f^ - ",) Xf - i)4 ^ i ^ - 2 ) 

(14) 

(Cont iniied) 


3 : 




(14) 

(Concluded) 





For the case of a two-bladed teetering rotor, the forces and moments may 
be obtained by setting e = 0, a^ = = 0 in the above equations. For a 

teetering rotor without cyclic pitch, as is used in many tail rotor systems, 
the forces and moments may be obtained by setting = 0. Further, 

since the tail rotor flapping frequency is much higher than that of the main 
rotor system, the tip-path plane dynamics may be neglected. Thus, for tail 
rotor applications, one simply sets a^ = aj = bj^ = 0, and ap = a^ - bj = 0 
in the above equations. The result is a set of basic quasi— static forces and 
moment expressions similar to those in the classical work (refs. 9, 10). 


These rotor forces and moments,, which were then transformed to the center 
of gravity of the aircraft, represent the contribution of the main rotor to the 
total aerodynamic forces and moments in the body axes of the aircraft. From 
figure 3 it can be seen that the transformation is simply: 


M, 


Tp sin ip - Hp cos ip 
-Tp cos ip - Hp sin ip 

L^p cos ip - Qp sin ip + Yphp + (Tp cos ip + Hp sin ip)dp 

Mj^P - (Tp 8ln ip •• Hp VOS lp)hp + (Tp cos ip + Hp sin lp)fp 

Qp cos ip + Lj^p sin I -i Ypfp + (Tp sin ip - Hp cos lp)dp 


(15) 
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Figure 3.- Main rotor forces and hub moments diagram with respect to 

center of gravity. 


aircraft 


^ ig'^’^^tion of the nine-degree-o f- freedom 
(DOF) mathematical model, which consists of three-DOF tip-path plane dynamics 

AY anrJr dynamics. Note that the aerodynamic forL terms Ax, 

AY and AZ, and the moment terms AL, AM, and AN are those contributed by the 
tail totor, fuselage, horizontal stabilizer and the vortical fin. 


CONTROL PHASING 


allv he diffr ? f <A2) that the flapping frequency will gener- 

nrl rhe r rotational frequency of the rotor system. There- 

exMM^ response to a cyclic-control input will no longer 

fho a phase. Proper control phasing will be required to achieve 

the desired flapping decoupling, that is, a longitudinal-control input 
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TABLE 1.- GENERALIZED SINGLE ROTOR HELICOPTER MATHEMATICAL MODEL 
(NINE DEGREES OF FREEDOM. FIRST HARMONIC) 


Dynamic equations: 


— r 

dt I . 


w = -I“^WjjI(i) + I~^M^(V,(»); ^,a; u) 


( 1 ) ■ C ::) 0 • 1. J 


Kinematic equations: 
Ve = [Lo 


Euler 

equations 


Tip-path- 

plane 

equations 


cos ({) 


-sin (p 


= 1 1 sin <t» tan 6 cos (p tan 0 


P sin <p sec 6 cos <p sec 0, 


where 


V = (u.v.w)’'’ ; 0 ) = (p.q.r)"^ ; a = (ao.aj.bj)*^ 

Vq = (Xg.ye.-h)^ 

T ^ 

F = mg(-sin 6, sin (|i cos 6, cos if) cos 0) ; u = (A, ,B, ,0„,9 ) 

— Ic lc O Ot 

F^ = (X,Y,Z)^ = (Xp + AX, Yp + AY, Z^ + AZ)^ 

M. = (L,M,N)^ = (Lp + AL, M„ + AM, N„ + AN)^ 


I = 0 


0 -I, 


0 -r q* 


= r 0 


cos 0 cos 


cos 0 sin i/) 


-sin 0 


^ 5 _£ sin 0 sin 4) sin — cos sin i|; sin 0 sin ({) sin ^ + cos 4> sin 4^ sin (p cos 0 
^In 0 cos (^cosil^-fslncpsin^^ sin0cos4>sin4/-sln(()cos4) cos p cos 0 



p 




produces only a steady-state longitudinal flapping, and a lateral-control input 
will produce only the lateral flapping response. A detailed discussion on the 
question of control phasing (or control mixing) is given in reference 8, 

A set of simplified equations based on hovering flight has been developed 
to achieve the desired flapping-response decoupling. The equations relate the 
required control phasing to the rotor system parameters, that is, 

3A, 


(16) 



35A 

36L 

9Bls 

'»ls 


36A 


8(5L 



3A, 


36A 

_ X VO/ 

1 - (4/3)e 

^Bls 

r-t 

1 

00 

1 

36A 

y[ 1 - (4/3)e] ‘^2 


_ 8(p2 - 1) . 

36L 

Y[1 - (4/3)e] 

3B, 

Is 

_ -fl- (8/3)e + 2 e2 


1 - (4/3)e 


07) 


with 


K eM„ yK 
p2 - 1 = P 




The so-called control advance angle Is given by 

3Bjg/36A 
3Aj g/96A 


tan 


P2 - 1 


y( 1/8 - t/3~+ e2/4) 


(18) 


The control crossfeed given by equation (17) has provided an 
for the entire low-speed flight regime. 


accuracy adequate 
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CONCLUDING REMARKS 




A simplified nine-degree-of-f reedom, main rotor, mathematical model has 
been developed for real-time, piloted simulation for exploratory Investiga- 
tions of helicopter flying qualities. The model explicitly contains the ' 

dynamic effects of the tip-path plane and several primary design parameters of ^ 

interest, such as flapping hinge restraint, effective hinge offset, blade Lock 
number, and pitch-flap coupling. Furthermore, being analytical, the model 
provides insight useful for designing simulation experiments. 

It should be emphasized, however, that the simplified mathematical model 
is suitable only for stability and control applications. Because of the basic 
assumptions used in its development, applications of the model should be 
limited to an advance ratio of about 0.3. Beyond that, the model must be 
modified to account for the stall and compressibility effects. Incorporating 
the rotor rotational degree of freedom, ground and turbulence effects, and 
inflow dynamics would enhance the model’s validity for a wider range of 
applications. 



APPENDIX A 


Flapping Equation of Motion 

Given the basic assumptions noted in the text, the flapping equations of 
motion have been developed in reference 8 for a rotor system containing design 
parameters, such as flapping hinge restraint, hinge offset, and pitch-flap 
coupling. For convenience, they are given here for both two-bladed teetering 
rotors and nonteetering rotors. 


Two-Bladed Teetering Rotor 
6 + -^ B + jj >2 + IH_ gin 2il> + — ^ — (1 - cos 

■ “'{[Kt «w - a 6w) <=»* i sin 6„) + I (£ 

+ a sin 6.,)Jsln * + ^20 cos + § sin B„) - sin B„ - i cos 6„) 
' 8 (fi t cos cos sin ’i')9o - | [(l + j p2)cos 4> 

- ~ cos 31/ijAj^ _ I jji + I yZjgin ,j, _ -L sin 


n COS 6, 


+ ^ (sin ^ (sin i(^)A 

where 


I 


(Al) 


.2 . . . S 


P-" = 1 + 


Ig«2 8 


N-Bladed Nonteetering Rotor 


®1 * 2*'[(J ■ 3 ‘ t) + “(j - s + £^)sln ♦jJbj 

" * I [<i - »1 - T (i - ' - #1" 

*1 (5 - = * t)<* ■ 2*pj| 

' “ *1 ■" ®”*n *'3 


(A2) 
(Continued) 
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■f (i- cos 6 „ + ^ sin S„)sin sin 6 „ - ^ cos 6 „)cos 

7 ? 1 <» - “<1 + P''> - l[(i " f) (l ■ ')“ *i 

4 (2 - ' * ■ ^[(l ■ ■5)'°® *1 * "(7 ■ I 

+ X (t • ^ * x)(cos (Tj - cos SI'jijAij - yj(j - §)sln *J 
+ y0 - f)(l - cos 2ip^) + ^ (j - £ + -y)( 3 sin 4;^^ - sin 34'^) 

+ t[(t ■ f) ■" "(1 ' I ♦i + ■¥ 0 ■ f)‘^ ■ 

* 2 [0 - 2) + “0 - ' + t)“‘" ♦i]*! + 10 ' ^)[( 

+ ^ Sin B„)sin 4 ^ ~ (n ^ cos S„^cos 4 J 

^)[(fi + fi 24^^ - sin 6 , 

- a cos 6„)sin 24 J 


li - 4llrg cos 6w 


where 


1 + 


°»6 , S , 

‘6 « 


i(i-^) 


i =■ 1, 2, 3, . . . , N 



(A2) 

(Concluded) 
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APPENDIX B 


Derivation of Rotor Force and Moment Expressions 

The derivations of the rotor forces and moments will be performed first 
in the wind-hub system (see fig. 2). Then, they will be transformed into the 
hub-body system to be added to the other components contributing to the total 
aircraft forces and moments. Only the development of thrust is given here to 
illustrate the procedure. The development of other forces and moments is 
given in a forthcoming report.^ 


The shear force for a single ith blade at azimuth is given by 

- mq + m(w - uq + pv) 


j^2pO^ + em^ + -^ qjcos 

+ em^ + pjsin 


where the aerodynamic force acting on the blade is again given by 

>R-e 


^a = 


I 


f (DR)2ac(U^20 + U^Up)dr' 

Substituting the following equations, 

- 

^'t OR “ e(l - cos 6) + y sin i|; + X cos 6 

S 


(Bl) 


(B2) 


(B3a) 


^P ^ ni “ ^ ^ ~ ^ ® tj; - ^ (x - e) 

+ *P cos S„ t a sin 6„)sli, ♦ + (- £ sin B„ + a J 


e = 0„ - A 


cos ^ sin + x6^ - KjB 


(B3b) 

(B3c) 


Into equation (B2) results in becoming 


2chen, R. T. N. ; and Decker, W. A.: Effects of Primary Rotor Parameters 

on Helicopter Stability and Control Characteristics (in preparation). 
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F. - I (aR)VRp + (1 - c)(l t J + „ sin ♦)„ sin *.j(6„ - cos t-Bi^sln*) 

■" [i ■"(I'"? 1 »)f Fl" »]et + |f (l-£2)+(i_c)n sin Jx 

~ ^ ~ ^ » sin ij()u sin + ||. (1 - ^2) 

+ (1 - e)M sin </-jp cos ’''je ~ f - f + f (l-e)2 sin B 
[i 2 (1- "^)sin cos 6„ + ^ sin B„)sln 

Q ®w + Q cos 6^) cos t|;j| ^ 

The flapping will be approximated by its tip-path plane representation: 


4 (1- e2) 


thus. 


8 - “ aj(t)cos \p - bj(t)sin ip 


e - Hq - (aj + bjQ)cos i|/ - (bj ~ ajfi)sin \(i 


(B5a) 


(B5b) 


6 - aj - (aj + 2bi« - 4> - (b^ - 2k^n - b^Q^)sln (B5c) 


Noting that 


" N 

Zj sin sin li-' = 0 

1=1 ^ i=l 1 

N N 

Jj cos t|» = cos il'! = 0 

i=l i=l i 


the thrust becomes 


^ ~ + mg - m(w - uq + pv) ] (] 

rnnrr 

obtained: following expression v 


18 


1 

(l/2)pacR(!T2R) 


N 


2 ■ "[t * T + 1 « - 


e^)X - 


^ 1^0 ®0 /I e 


i=l 


(i- 


where 


+ fjy^d - e) 0 (, + -^ (1 - e 2 )q^ _ ^ 

- [a(,p 2 (i _ e) _ b^y(l - e 2 )jk^ + ^ . 

+ ^ (bi - ajfl) I (1 - e )2 

+ f (1 - e 2 )^g cos 6 „ + ^ sin + 0 (N) 


0 ( 4 ) = 0 


(for four-bladed rotor) 


0(3) ~ J [m2(i _ e)Bj^ - bjy2(i _ + a^u^Cl - e)]sin 3ip 


2 

+ 4 [U^(l - e)Aj^ - (ajKj + bj)y2(i . ^)|cos 3«|; 


(for three-bladed rotor) 


Finally, the thrust may be expressed as 

T ^ N 


- I PacK(.R) 2 ji (1 . . [i . ^ (X - ,)]e„ . [i . ^ (1 - e^)]e. 

(1 - - a„ + L (1 . j"! 

*b. [i (1 - j| 

- =')(§ =»* 8 v + t am 6 „) * 




+ nig - m(w - uq 


where 


0 ( 4 ) = 0 


(for N = 4 ) 


0(3) ^ [i,2(i _ _ bj^2Q _ ^ a^u^(l - e)]sin 34> 

2 

+ 4 [u^a - e)Aj^ - u2q _ e)(a^K, f bj)]cos 3.^ 


(for N = 3) 
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1 ¥ 


)] 

>»ic ; 

e 2 ) 

r 

( 87 ) 


pv)] 


( 88 ) 


By dropping the high harmonic contributions, the thrust can now be expressed as 


T = I pacR(nR )2 I (1 - e2)^ [3 ^ " e)J Og + 


+ V 


- ^ (1 - - Kjbj) - at 


j^3 ''■2 ^(1 - e)j 


- N[agMg + mg - m(w - uq + pv)] 


cos ^ sin 3^, 


P'^^’ generally small compared 

o other terms, they are dropped, as was shown earlier in equation (2). 
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lb Abstract 

This report documents a simplified analytical mathematical model of 
the helicopter main rotor; the model was developed primarily for real-time 
pilot— in— the-loop investigation ot helicopter flying qualities. The mathe- 
matical model explicitly Includes the tip-path plane dynamics and several 
primary rotor design parameters, such as flapping hinge restraint, flapping 
hinge offset, blade Lock number, and pitch-flap coupling. The model has 
been used in several exploratory studiea, recently performed at Ames 
Research Center, of the flying qualities of helicopters with a variety of 
rotor systems. 

The report describes the basic assumptions used and the major steps 
involved in the development of the set of equations listed. The equations 
consist of the tip-path plane dynamic equation, the equations for the main 
rotor forces and moments, and the equation for control phasing requited to 
achieve decoupling in pitch and roll due to cyclic inputs. 
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